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ABSTRACT 



According to a recent modified model of gravity at large distances, a radial 
constant and uniform extra-acceleration Aft in = of Rindler type acts upon 

a test particle p in the static field of a central mass M if certain conditions are 
satisfied. Among other things, it was proposed as a potentially viable explanation 
of a part of the Pioneer anomaly. We study the impact that an anomalous 
Rindler-type term as large as |ARj n | ~ 10~ 10 m s~ 2 may have on the the orbital 
dynamics of a typical object of the Oort cloud whose self-energy is quite smaller 
than its putative Rindler energy. By taking a typical comet moving along a highly 
eccentric and inclined orbit throughout the expected entire extension of the Oort 
cloud (~ 0.02 pc — 1 pc), it turns out that the addition of an outward Rindler-like 
acceleration, i.e. for A Rin > 0, does not allow bound orbits. Instead, if A Rin < 0, 
the resulting numerically integrated trajectory is limited in space, but it radically 
differs from the standard Keplerian ellipse. In particular, the heliocentric distance 
of the comet gets markedly reduced and experiences high frequency oscillations, 
its speed is increased, and the overall pattern of the trajectory is quite isotropic. 
As a consequence, the standard picture of the Oort cloud is radically altered 
since its modified orbits are much less sensitive to the disturbing actions of the 
Galactic tide and nearby passing stars whose effects, in the standard scenario, 
are responsible for the phenomenology on which our confidence in the existence 
of the cloud itself is based. The present analysis may be supplemented in future 
by further statistical Monte Carlo-type investigations by randomly varying the 
initial conditions of the comets. 

Subject headings: Experimental studies of gravity; Experimental tests of gravitational 
theories; Modified theories of gravity; Oort cloud; Orbital and rotational dynamics 
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PACS: 04.80.-y; 04.80.Cc; 04.50.Kd; 96.50.Hp; 96.25.De 



1. Introduction 



Recently, iGrumillerl (120101 ) . under certain assumptions, put forth a quite general model 
for the gravitational field of a static central object of mass M at large distances r from 
it. As a result, the acceleration felt by a test particle p in the field of M turns out to be 
modified by the appearance of a Rindler-type additional term. Thus, the total acceleration 
becomes (IGrumillerl |2010| ; ICarloni et al.l 1201 ll ; iGrumiller fc Preisl 1201 ll ) 



GM 



A 



Rin; 



where G is the Newtonian gravitational constant, and 



Rin 



.4 



(2) 



the sign of Ajt\„ is left indeterminate by the theory. Notice that IGrumiller fc Preisl ( 120111 ) 
and ICuletuI (I201lh assumed that A^ in is a universal co nstant. Concerning its sign and 
magnitude, IGrumiller Sz Preisl (120111 ) and ICuletuI (120 111 ) took 



\A 



Rin 



1 x 10" 10 m s 



(3) 



Grumiller &: Preisl (120111 ) assumed an inward direction. 



Importantly, ICarloni et al.l (120111 ); IGrumiller fc Preisl ( 1201 ll ) argued that, if eq. (JTJ) has 
to be valid for a given physical system, then the following condition 



Gm 



(I 



V -<\A 



Rin | T 



(4) 



must be fulfilled. In eq. fll]) m p and d p are the mass and a typical size of the test 
particle p, respectively. Otherwise, the self-energy of the test particle would overwhelm the 
Rindler energy, and the consequent particle's backreaction on the background would not 
be negligible. In this case, it would not be possible to assume the universal (=maximal) 
value for the Rindle r acceleration. For more details, see the discussion in Section V of 
Carloni et al.l ( 120111 ). Actually, it is easy to show that eq. is not valid for the planets of 
the solar system since it is 



Gm p 
dp \A Kin \r 



(6 x 10 2 - 1 x 10 5 ) kau 



(5) 



where 1 kau = 1000 ast ronomical units = 4.85 x 10 3 pc. ICarloni et al.l (120 111 ); 
Grumiller fc Preisl ( 1201 ll ) noticed that eq. (jlj) is, instead, fully satisfied for the Pioneer 
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10/1 1 probes, affected by the well known Pioneer anomaly ([Anderson et al.l 11998 



200 21) which may partly be explained by th e Rindler-type acceler ation (IGrumiller 



2010 



Carloni et al.l 12011c IGrumiller fc Preisl l201lh . ICarloni et al.l ( 1201 ll ) obtained the following 
constraint 

\Amn\ < 3 x 10~ 9 m s~ 2 (6) 

from the propagation of electromagnetic waves, for which the previous caveat is not a 
concern. 

Actually, the solar system host several n atural objects orbiting the Sun for which 
eq. (|4]) holds: they reside in the Oort cloud (IOpik!ll932l ; IOortlll950l ). It is a reservoir 
of frozen cometary nuclei which is supposedly located in the remote peripheries of the 
solar system. It should be a remnant of the early stages of the formation of the solar 
system, and it likely formed as a consequence of scattering of planetesimals by the giant 
planets (iDuncan et al.l 19871 ; iHiguchi et al.ll2006l ). The Oort cloud has likely a spheroidal 
shape (jWeissmarJ 119961) a nd a size ranging from about 5 — 10 kau up to 150 — 200 kau 
(ILevison Sz Donnesl 120071 ). Its existence was conjectured by no ticing that the typical 
lifetimes of comets near the Sun is of the order of 10 4 — 10 5 yr ( Levison fc Duncan 19941 ) 
due to non-gravitational (sublimation) and gravitational (strong interactions with planets) 
phenomena, while the age of the solar system is 4 — 5 orders of magnitude larger. Thus, a 
continuous resupplying from a much more remote, longer-lived source should take place in 
order to maintain the cometary population in a steady state, as it is observed. Long-perio cjj] 



(Ph > 200 yr) comets, characterized by highly i nclined and eccentric orbits (Morbidellil 



20051 ). would originate just from the Oort cloud (jMorbidellil 120051 ; lDuncanll2008l ). They are 
injected into observable orbits in the p lanetary reg ions of the solar system by g ravitation al 
inter a ctions with interstellar me dium (ISternlll990l ) and giant molecular clouds (jMazeeva 
20041 ; Ijakubfk fc Neslusanl l2008t l , near by passing stars (lHillslll98ll: lBobvlevll2010al lbh and 
the Galactic tide (radial an d vertical') (|Rickman et al.ll2008t iMasi et al.ll2009l ): as a result, 



comet showers may occur (IHeisler et al.lll987l ; iMatese et al.lll995l ). Also other existing 



minor objects of the solar system like Centaurs, highly elliptical tran s-Neptunian objects 



and J upiter-family comet population may come from the Oort cloud ( lEmelyanenko et al. 



20071 ). Connections of terrestrial cratering with such co metary showers originating from the 
Oort cloud have been investigated by some researchers (jWickramasinghe fc Napierll2008l ). 

By taking the Hallej|§ c omet as representati ve, a typical Oort object p may be thought 



as characterized by a mass (jCevolani et al.l 119871 ) 



m r 



2.2 x 10 14 kg 



(7) 



1 The thresho ld of 200 yr was chosen mainly for historical reasons: it is arbitrary 
jMorbidellil l2005| ). 

2 Although Pb = 75.3 yr, its orbital c haracteristics po int towards a capture from the 
long-period population of the Oort cloud ( IFernandezI 120021 ) . 
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and mean size (iLamy et al.l 12004 ) 



dp = 11 km, 



with a mean density of (ISagdeev et al.lll988l ) 

p p = 0.6 g cm~ 3 , 



(9) 



although uncertainties in several parameters and assumptions (jPeald Il989f ) may lead to a 
range of values as large as 



p p = 0.2 — 1.5 g cm" 



Such figures and eq. ([3]) imply that 



0.0892 au 



d p \A Rin \ r 



(10) 



Thus, the condition eq. (j4|) is fully satisfied over the entire extension of the Oort cloud also 
for bodies which can be much larger and denser than a typical cometary nucleus. Indeed, 
eq. (jlj) can approximately be posed as 



nGppdl 

I A I ~ ' 



(12) 



for, s ay, p p = 5 g cm 3 , which is the mean densitjH of the rocky planets ( McFadden et al. 



20071 ). and d p = 250 km we would have 



HP p - = 4.4 kau. 



A 



Rin| 



(13) 



In this paper we will explore the consequences that the existence of a Rindler-type 
extra-acceleration, with the characteristic of eq. (TS]) and eq. (EJ), would have on the orbital 
motion of a typical Oort cloud object. Should the resulting orbital pattern radically 
differ from the standard Newtonian one, shadows on the Rindler-like acceleration would 
be casted since the entire dynamical history of the Oort cloud should be re-written and 
all the inferences nowadays accepted, based on the Newtonian picture of the Oort cloud 
and of its interaction with the surrounding stellar and Galactic environment, would not 
be valid anymore. Actually, non-negligible modifications of the Newtonian orbits are 
expected since the mean Newtonian accelerations of an Oort comet may range from 
v4 N = 2 x 10~ 10 m s~ 2 (r = 5 kau) to A]<s — 1.5 x 10~ 13 m s~ 2 (r = 200 kau). Concerning 
general relativity, it will not be considered in the rest of the paper since its effects are 



3 Clearly, it is highly unlikely that such dense bodies can really exist in the Oort cloud. 
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totally negligible in the present context. Suffices it to say that the 1PN term of order 
0(c~ 2 ) causes a comet acceleration as small as 



A 1PN ~ < 1(T 21 m s- 2 , (14) 



where c is the speed of light in vacuum. As we will see in Figure [21 (Section 12. ip . no special 
relativistic effects will come into play since, although generally increased by the Rindler 
acceleration, the speed of the Oort comet will remain several orders of magnitude smaller 
than c. 

The plan of the paper is as follows. In Section [2] we numerically compute the trajectory 
of a typical Oort object for inward (Section 12.11) and outward (Section I2.2p directions 
of Ap{j n . In Section [3] we investigate if an Oort comet acted upon by a Rindler-type 
acceleration is sensitive to the perturbing actions of the Galactic tide (Section 13. ip and of 
nearby passing stars (Section 13. 2p which, in the standard Newtonian scenario, cause the 
phenomenology on which our confidence in the existence of the Oort cloud is based. In 
Section 0] we summarize and discuss our findings. 



2. Numerically produced orbits of Oort comets affected by the Rindler 

acceleration 

2.1. An inward Rindler acceleration 

The orbital effects of eq. (l2l)-eq. ([3]) on an Oort comet p cannot be computed 
perturbatively since, according to eq. (j3J), the magnitude of Am Q is of the same order of 
magnitude of, or larger than the Newtonian acceleration experienced by p. 

Thus, we numerically integrate the modified equations of motion of p in a heliocentric 
frame endowed with cartesian coordinates. We adopt initial conditions, shown in Table 
[H such that, in the limit Ani n —> 0, the resulting orbit reduces to a standard Keplerian 
ellipse. It is characterized by large values of its eccentricity e and inclination / to the 
reference {x,y} plane, and the semi-major axis a is chosen in such a way that the comet's 
motion covers almost all the expected extension of the Oort cloud. The other standard 
Keplerian orbital elements are the longitude of the ascending node i?, the argument of 
pericenter u and the true anomaly /. In order to facilitate a comparison with other studies 
on the Oort cloud, we adopted as reference frame a heliocentric inertial one having the 
x axis directed towards the Galactic Center, the y axis direc ted along the Sun's G alactic 
velocity, and the z axis completing a right-handed system@ ( Fouchard et al. 2005 ). The 



4 The Galactic plane is tilted by 62.6 deg to the celestial equator (jSullivan!ll984l ). so that 
it is almost at right angle to the ecliptic. 
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Table 1: Initial conditions, in kau and kau Myr -1 = 4.74 m s _1 , adopted for the numerical 
integration of the modified equations of motion of an Oort comet acted upon by the Rindler- 
type acceleration of eq. ([2]). In terms of standard Newtonian mechanics, they correspond to 
a Keplerian ellipse with a = 102.507 kau, e = 0.955, I = 77 .41 deg, Q = -139.87 deg, uj = 
321.94 deg, / = —152.14 deg. They c an be thought as referred to the heliocentric fixed frame 
of Figure 1 in 
distance is q = 
kau. 



Fouchard et al.l (120051 ) . The orbital period = 32.81 Myr. The perihelion 
a(l — e) = 4.606 kau, and the aphelion distance is Q = a(l + e) = 200.407 



Xq (kau) |/o (kau) z (kau) Xq (kau Myr r ) y (kau Myr : ) i (kau Myr ^ 



45 



35 



10 



-23 



■15 



-15 



results of the integrations are displayed in Figure [TJ and Figure |2] in which the features of 
the Keplerian ellipse corresponding to the initial conditions of Table [TJ are shown in blue 
for comparison. The qualitative differences with the Newtonian case are striking: bound 




Fig. 1. — Sections in the coordinates planes of the numerically integrated orbits for the 
Oort comet of Table [TJ over one Keplerian orbital period P\,. The dash-dotted red lines are 
the projections of the modified orbit by assuming that the Rindler acceleration is directed 
towards the Sun. The dashed blue lines are the projections of the Keplerian ellipse. 

trajectories still occur, but they radically differ from the Newtonian ones. It is important 
to note that the condition of eq. (Hj) is always fulfilled, as shown by the right panel of 
Figure [2 This is a relevant point since, in principle, it may happen that the comet enters a 
region in which its gravitational self-energy becomes larger than its putative Rindler energy, 
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Fig. 2. — Heliocentric distances, in kau, velocities, in units of c, and ratio Gm p / (d p \A-m n \ r) 
of the numerically integrated orbits for the Oort comet of Tabled] over one Keplerian orbital 
period P^. The dash-dotted red lines refer to the modified orbit by assuming that the Rindler 
acceleration is directed towards the Sun. The dashed blue lines refer to the Keplerian ellipse. 



thus destroying the validity of the approach followed which was implemented by keeping 
^Rin 7^ throughout the whole integration. As anticipated in Section [T] the non-relativistic 
treatment is justified by the smallness of the comet's speed v, displayed in the middle panel 
of Figure [2j Globally, the heliocentric distance r is strongly reduced, especially in the 
regions which in the Newtonian case correspond to the aphelion: the minimum distance 
suffers a relatively smaller reduction with respect to the Newtonian case. Indeed, the left 
panel of Figure [2] tells us that, in the Rindler case, r oscillates with high frequency between 
about 2 kau and 50 kau. 



2.2. An outward Rindler acceleration 

If the Rindler-type acceleration is radially directed from the Sun to the comet, no 
bounds orbit may exist, as depicted by Figure |3] which refers to the same initial conditions 
of Table [TJ We adopted decreasing values of A Rin with respect to |A Rin | = 1 x 10~ 10 m s~ 2 , 
but the result is substantially the same. 



3. Consequences of a Rindler-type acceleration on a perturbed Oort cloud 

It seems plausible to expect that the Rindler trajectories are less sensitive than the 
Newtonian one to those disturbances which affect the Oort cloud yielding those phenomena 
which are the basis of our confidence in the existence of the cloud with the features usually 
attributed to it. If so, the existence of the Rindler acceleration would drastically alter the 
dynamical history of the Oort cloud. 
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Fig. 3. — Numerically integrated heliocentric distances r, in kau, of the modified trajectory 
(dash-dotted red lines) of an Oort comet with the initial conditions of Tabled] for decreasing 
positive values of A Rin .The duration of the integration is one Keplerian orbital period P^. 

3.1. The Galactic tide 



The effect of the Galactic tide on the orbit of an Oort comet ( jHeisler fc Tremaine 



19861) can be obtaine d from the following expression of the Galactic tidal acceleration 
flFouchard et~alll2005h 

4 (tid) _ 



A 



(tid) 



.4 



(tid) 



—Q\x cos (Cl Q t) + Q 2 y sin (Cl t) 
—Q\x sin (Clot) — G 2 y' cos (Clot) 



(15) 



The various quantities entering eq. (fl5|) are 

x = x cos (Clot) + V sin (Cl t) , 

y = —x sin (Cl t) + y cos (Cl t) 
Gi = -(A-B)(3A + B), 

g 2 = (a — b) 2 , 

Gs = 4:irGpo - 2 (B 2 - A 2 ) , 



(16) 



and 



(17) 
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where A and B are the lOortl (119271 ) constants, po is the mass density in the solar 
neighborhood, and flo is the frequency of solar revolution^] around the Galaxy. Both 
visible stars and Dark Matter concur to the local mass density in the solar neighborhood. 
Co ncerning the bar yonic stellar component, the local mass density of main-sequence stars 
is flReid et al.lliooih 

Pms ~ O.O31M pc" 3 = 2.09 x 10" 24 g cm" 3 , (18) 

while, according to th e most recent results, th e Dark Matter local density in the Galactic 
solar neighborhood is ( Ide Boer fc Weberll201l[ ) 

p DM = 1.3 Gev cm" 3 = 2.31 x 10" 24 g cm" 3 . (19) 

Thus, we will assume 

po ~ 4 x 10" 24 g cm" 3 = 1.5 x 10 22 kg kau" 3 . (20) 



The frequency of t he Galactic revol ution of the Sun Q can be evaluated from its 
circular rotation speed (IReid et al.l 120091 ) 0o = 254 km s _1 and its distance to the Galactic 
center flReid et al.ll2009h R = 8.4 kpc as 



\n \ = % = 0.0309 Myr" 1 . 
R 



Thus, its period of revolution is 



T = 203.176 Myr. 



Concerning the Port constants A and B entering eq. ([T 
accurate to 5%, are (IFeast fc Whitelock!ll997l ) 

A = 14.82 km s" 1 kpc" 1 = 0.0151 Myr" 1 



(21) 
(22) 

their most recent values, 



B 



■12.37 km s" 1 kpc" 1 = -0.0126 Myr" 1 . 



(23) 



Thus, eq. (JUD and eq. 



inserted in eq. (H7J), yiekf] 
g 1 = -9.1261 x 10" 4 Myr" 2 , 



Q 2 = 7.7325 x 10" 4 Myr" 2 , 
Q 3 = 3.9626 x 10" 3 Myr" 2 . 



(24) 



5 Since the motion of the Sun about the Galactic Center is clockwise in the frame adopted 
( IFouchard et al.l 120051 ). f2 is negative in it, i.e., it is directed along the — z axis. 

6 Cfr. with the values by iLevison et al.l (|200l[ ). 
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In Figure S] we show the numerically integrated Rindlerian and Newtonian trajectories, 
including eq. (I15p . over a full revolution of the Sun around the Galaxy. We used the same 
initial conditions of Tabled] for the Oort comet. As expected, the differences between the 



A Rin = -1 10- 1D ms- 2 




y (kail) I (Myrl 

Fig. 4. — Coordinate sections of the numerically integrated trajectories and heliocentric 
distances r of an Oort comet affected by the Galactic tide of eq. ffl5l) with (dash-dotted red 
lines) and without (dashed blue lines) the Rindler-type acceleration of eq. 02])- The initial 
conditions of the Oort comet are those of Table [TJ The duration of the integration is one 
full Galactic revolution (T = 203.176 Myr). 

Newtonian and the modified scenarios are remarkable: a comparison with Figure [1] and 
Figure |2] shows that the Rindler trajectory is left substantially unaffected by the Galactic 
tide, contrary to the Newtonian one. It can be shown that also for an outward Rindler 
acceleration the Galactic tide does not have influence in the sense that the comet does not 
remain bound, as in Figure El 



3.2. A close encounter with a passing star 

Here we consider the perturbing action of a passing star s with m s = M over an 
unperturbed Keplerian orbital period Pb of the Oort comet. The initial conditions of s are 



listed in Table [2j In the numerical integration of the equations of motion we assume that 



Table 2: Initial conditions, in kau and kau Myr -1 = 4.74 m s , of a perturbing passing star 
s adopted for the numerical integration of the modified equations of motion of an Oort comet 
acted upon by the Rindler-type acceleration of eq. ([2]). They correspond to a heliocentric 
distance = 3.2 pc and a speed u„ = 0.61 km s -1 . 



Xq (kau) (kau) 


4 S) (kau) 4 S) ( kau Myr- 1 ) 


yj s) (kau Myr" 1 ) 


4 s) (kau Myr" 1 ) 


270 370 


470 -42.18 


-63.28 


-105.47 



s moves uniformly with v®(t) = v$\ The result is depicted in Figure El It can be noticed 
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Fig. 5. — Coordinate sections of the numerically integrated trajectories and heliocentric 
distances r of an Oort comet affected by the passing star of Table [2] with (dash-dotted 
red lines) and without (dashed blue lines) the Rindler-type acceleration of eq. ([3]). The 
initial conditions of the comet are those of Table [TJ The duration of the integration is one 
(unperturbed) Keplerian period Pb- 

that, while the Keplerian ellipse is, actually, distorted by the Newtonian interaction with 



'In this particular case, the comet is stripped away by the stellar passage: indeed, this is 
one of the causes of the gradual depletion of the Oort cloud (IRemy fc Mignardlll985l ). 



13 



s, the trajectory computed with the Rindler acceleration remains substantially unaffected. 
Also in this case, it can be shown that no bound orbits can exist if Ani n is directed outward. 



4. Discussion and conclusions 

The analysis performed in this paper should be considered just as a preliminary one. 
It aims at exploring semi-qualitatively how the presence of an additional radial Rindler-like 
acceleration may affect the orbital motion of bodies whose self-energy is smaller than 
their putative Rindler energy like the objects moving in the Oort cloud. We only choose 
a specific set of initial conditions which correspond, in the Newtonian scenario, to a 
typical orbital configuration for an Oort comet moving on a highly eccentric and inclined 
ellipse which almost extends throughout the entire expected extension of the Oort cloud. 
More refined numerical analyses may explore larger ensembles of initial conditions with 
a statistical approach. It would also be interesting to repeat Monte Carlo simulations of 
th e dynamical evo l ution of th e Oort cloud over the age of the solar system like, e.g., those 



by iRickman et al.l (120081 1 and iMasi et al.l (120091 ) by explicitly including a Rindler-type 



extra-acceleration, and inspect how several key features of the cloud change with respect to 
the usual Newtonian scenario. 

In our numerical analysis we found that the standard Newtonian picture is notably 
altered by an additional radial Rindler acceleration as large as 10 -10 m s -2 which, in 
the case of an Oort comet, would not be smaller than the Newtonian one. Bound orbits 
are not possible if AR in is directed outward, even for Ain Q = 10~ 13 m s~ 2 . Conversely, if 
the Rindler acceleration is directed from the comet to the Sun the resulting trajectory 
is limited in space, but it is completely different from a Keplerian ellipse. In particular, 
the heliocentric distance r is greatly reduced, and it experiences high-frequency variations 
during a Keplerian orbital period P^. Moreover, the spatial pattern of the modified 
trajectory is quite isotropic over P^. As a consequence, the Rindlerian trajectory is much 
less sensitive than the Newtonian one to disturbing effects like the Galactic tide and nearby 
passing stars. 

Although, as pointed out above, more extended numerical investigations should be 
implemented by varying the initial conditions in a Monte Carlo fashion, it is difficult to 
believe that the main features of the Oort cloud, which lead to a general consensus about 
its existence, may be preserved by the existence of a Rindler-type extra-acceleration of 
the order of 10~ 10 m s -2 . Our results are not necessarily limited to the specific model by 
Grumiller indeed, they are quite general and are valid for any hypothetical constant 



and uniform radial acceleration. In the case of the Pioneer anomaly, the magnitude of the 
anomalous acceleration may be as large as |Api Q | < 1 x 10 -9 m s -2 and all the previous 
results remain qualitatively valid, being the anomalous behavior of an Oort comet even 
more remarkable from a quantitative point of view. Finally, we notice that using the Oort 
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cloud may, in principle, be useful also for other long-range modified models of gravity. 
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